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Lunar gamma ray emission seen during the first year by Fermi 
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N.Giglietto on behalf of the Fermi Large Area Telescope Collaboration 
INFN Bari and Dipartimento Interateneo di Fisica- Universitd e Pohtecmco di Bari 

We report the detection of the lunar gamma-ray emission during the first year of Fermi-LAT 
observations. Such emission is produced by cosmic ray nuclei interacting with the lunar surface. 
Thanks to the solar minimum conditions and the reduced effects of heliospheric modulation, the 
lunar flux was at its maximum due to the increased flux of Galactic cosmic rays hitting the lunar 
surface. Fermi-LAT instrument has a superior sensitivity, angular resolution, and observes the whole 
sky every two orbits, ft is the only gamma-ray mission capable of detecting the lunar emission with 
high confidence and to monitor it over the full 24th solar cycle. We also report the status of a search 
of the gamma-ray emission from major planets and asteroid populations in the ecliptic plane. 



I. INTRODUCTION 

The 7-ray emission produced by solid solar system 
bodies is due to the interactions of Galactic cosmic ray 
nuclei (mainly protons) with their surface layers. The 
main processes involved are the production and de- 
cay of TT^s and kaons by ions, bremsstrahlung by elec- 
trons and Compton scattering of the secondary pho- 
tons. The 7-ray telescope EGRET on the Compton 
Gamma-Ray Observatory (CGRO), operated from 
1991 to 2000 and detected the 7-ray emission from 
the Earth the Moon [3, [1], and the Sun 0. Nu- 
clear 7-ray emission produced in the lunar surface in 
low-energy cosmic ray and secondary neutron reac- 
tions was observed by the Lunar Prospector [4] and 
used in the analysis of the composition of the regolith. 
Although similar physical processes are involved, the 
7-ray spectra of the Earth, the Moon, and the Sun are 
very different. The Moon is so far the only observed 7- 
ray emitting body with the solid surface. For the Sun, 
the 7-ray emission from the disk, due to the interac- 
tions of cosmic ray nuclei with the solar atmosphere 
d, is accompanied by extended and brighter 7- 
ray emission due to the inverse Compton scattering of 
Galactic cosmic ray electrons off solar photons [2, 6, 7] . 

Early analysis of EGRET observations of the Moon 
yielded the integral flux of F{E > 100 MeV) = 
(4.7 ± 0.7 syst) x lO"'^ cm^^ 3. A later re- 
analysis confirmed the detection and yielded a flux 
F{E > 100 MeV) = (5.55 ± 0.65) x 10"^ cm'^ 
averaged over the entire mission duration Q . 

Calculations of interactions of cosmic rays with the 
lunar surface are fairly straightforward and involve a 
well-measured spectrum and composition of cosmic 
rays near the Earth and composition of the Moon 
rock, the regolith, which was quite rigorously stud- 
ied using the samples returned by the lunar missions 
as well as by the remote sensing. The first calcula- 
tion of the lunar 7-ray emission was done by Morris Q 
using the cross section data and techniques available 
at that time. Recent detailed flux calculations^ de- 
pends on the level of heliospheric modulation and is 
generally consistent with observations and were done 
using the Geant4 framework [l3|. It was shown that 



the spectrum of 7-rays from the Moon is steep with 
an effective cutoff around 3-4 GeV (600 MeV for the 
inner part of the lunar disk). Due to the kinematics 
of the collision, the secondary particle cascade from 
cosmic ray particles hitting the lunar surface at small 
zenith angles develops deep into the rock making it 
difficult for 7-rays to get out. Therefore the lunar 7- 
ray emission is produced by a small fraction of splash 
particles in the surface layer of the moon rock. High 
energy 7-rays can be produced by cosmic ray parti- 
cles hitting the Moon surface with a more tangential 
trajectory; thus only a very thin limb contributes to 
the high energy emission. 

A similar emission should be detected from any 
other solid object in the solar system. Therefore 
planets should emit 7-rays produced from pion-decays 
coming from the hadronic interactions by cosmic-rays 
hitting the surface or the atmosphere of these bod- 
ies. Recent papers [ll|, [l^ have studied interactions of 
cosmic rays with populations of small bodies in the 
solar system. In particular, gamma-ray emission from 
asteroids, rocks, and dust at the outskirts of the so- 
lar system, e.g., the Oort Cloudp^, may contribute to 
the isotropic gamma-ray b ackg round. The gamma-ray 
emission from the Moon0, [ij] can serve as a template 
for such studies. 

We report here the updated observations of the lu- 
nar 7-ray emission, previously presented [l^ [3, 
and the status of a search for emission from smallest 
solid objects, such as Jovian troians||19i]. 



II. DATA SELECTION 

The data used in this study were collected during 
the first year beginning from 4 August 2008 during the 
solar minimum of activity. We have applied a zenith 
cut of 105° to eliminate photons from the Earth's 
limb, and used the "Diffuse" class p^, correspond- 
ing to events with the highest probability of being 
true photons. The analysis is performed using the 
standard maximum-likelihood spectral estimator pro- 
vided by the LAT Science Tools [25| package (version 
v9rl5), which is available from and using p6_v3 post- 
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FIG. 1: Count map in a Moon centered frame, Right Ascension and Declination offsets in degrees respect to the Moon 
position in abscissa and ordinate respectively. The offset ranges are ± 8° for Right Ascension and ± 5° for Declination. 
The image has been obtained by using photons with E> 100 MeV, a bin width of 0.2 degrees with a gaussian smoothing 
2 bin radius. The colour scale is linear with the counts. 



launch Instrumental Response Functions (IRFs). We 
also used the standard Science Tools provided by the 
Fermi Science Support Center, version v9rl5 and IRFs 
(Instrumental Response Functions) version P6_V3. 

The position of the Moon has been computed using 
a JPL[2J] library interface. Photon events for analysis 
were selected in a frame, 15 degree radius, and cen- 
tered on the istantaneous lunar position. We used the 
unbinned likelihood analysis technique typically used 
for astrophysical sources but, since the Moon moves 
quickly across the sky, we had to take additional pre- 
cautions to determine the emission spectrum. In or- 
der to avoid strong variations of background photons, 
we excluded time intervals when the Moon was close 
to the galactic plane or bright sources. We there- 
fore required that the Moon was at least 30° from 
the galactic plane i.e. \Bmoon\ > 30°, and we remove 
also any time intervals in which any individual bright 
object has an angular distance less than 5° from the 
Moon. The real difficulty in obtaining a consistent 
lunar 7-ray flux is the variation of the diffuse galactic 
background and extragalactic emission as the Moon 
moves through the sky. Therefore to evaluate the dif- 
fuse background in proper way we use the fake source 
method [23. [23j by defining a "fake" moon that fol- 
lowed the lunar trajectory but 30 degrees displaced 
along the true trajectory. The "fake" Moon was there- 
fore exposed to the same celestial sources as the true 



Moon and the event observed in the frame centered 
on the fake Moon make a good description of the dif- 
fuse background. Moreover we verify that this back- 
ground estimation is consistent with the background 
flux values measured using larger values of angular 
displacement along the lunar trajectory to define the 
fake source. 



III. LUNAR OBSERVATION RESULTS 

Fig. [1] shows the count map of photon events with 
E> 100 MeV in offsets of celestial coordinates relative 
to the Lunar position. Emission from the Moon is 
clearly visible and centered on the expected location 
in this relative coordinate frame. 

Fig. [5] shows the count map of the events with 
E> 100 MeV and within a 15° from the Moon cen- 
ter projected onto Right Ascension while in Figl3]the 
events are projected onto Declination. The coordi- 
nates are offsets in degrees of celestial coordinates rel- 
ative to the Moon position. In these figures the counts 
observed for the "fake" Moon are superimposed and 
demonstrate the need to carefully consider the back- 
ground before any lunar analysis can be performed. 
Simulated data confirm that the observed shape is 
in agreement with a pointlike source and the calcu- 
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FIG. 2: Count map of photons having an angular distance 
less than 15° from the Moon center and with E> 100 MeV, 
as a function of Right Ascension offsets in degrees respect 
to the lunar nominal position. Superimposed as dashed 
line the fake Moon count map distribution. Observed 
data are consistent with the expected angular resolution 
for 100 MeV photons. 



FIG. 3: Count map of photons having an angular distance 
less than 15° from the Moon center and with E> 100 MeV, 
as a function of Declination offsets in degrees respect to the 
lunar nominal position. Superimposed as dashed line the 
fake Moon count map distribution. Observed data are con- 
sistent with the expected angular resolution for 100 MeV 
photons. 



lated point spread funtions[26| of the LAT[20] for the 
expected lunar gamma-ray spectrum. In a more de- 
tailed analysis [215 of a one- year database we will per- 
form carefully the analysis of the shape of emission to 
determine its true extent and whether we can discern 
evidence for the expected limb brightening. 

To study the spectrum of the gamma-ray emis- 
sion we have used the standard unbinned maximum- 
likelihood spectral estimator provided with the LAT 
science tools gtlike. This preliminary analysis was 
performed by fitting the "fake" moon data to model 
the background and the Moon data sample with ei- 
ther a simple power law or other functional forms like 
a log-parabola summed to the background model. 

The fitted values obtained using a simple power law 
give a good test-statistics value of 7320.3 (defined as 
TS = 2(logL — logLo); being L and Lq the likeli- 
hood values when the source is considered or not) and 
indicate a power law index of —3.13 ± 0.03. As a 
result of these fits we estimate the observed fiux as 
F{E > 100 MeV) = 1.1 ± 0.2 x lO'^ cm-^s"!. The 
error includes the estimation of overall systematic er- 
ror of about 20% for these measurements, due essen- 
tially to the uncertainties of the detector effective area 
and due to the inefficiencies in gamma-ray detection 
due to pile-up effects from near coincidences with cos- 



mic rays in the LAT detector. A comparison of our 
data with models [i| are in good agreement with the 
shape of the expected spectrum and compatible with 
the levels computed of solar activity. The detailed 
analysis of the differential spectrum of the lunar emis- 
sion, in particular below 100 MeVfol. [l3|. is currently 
in progress. 

Our measurement of the lunar 7-ray flux is about 
twice the average fiux measured 5.5 x 10~^ cm~^ s"""^ 
previously obtained by EGRET [^]. However we know 
that the EGRET data subdivided into two subsam- 
ples, data taken during the solar maximum 1991 - 
1992 and those taken during period of moderate ac- 
tivity 1993 - 1994, yield to 3.5 x 10"'^ cm-^ s"! and 
7 X 10~^ cm~^ s~^ respectively, demonstrating that 
the flux is higher near the solar minimum. Therefore 
the Fermi-LAT observations reported in this paper, 
that are taken entirely during the period of solar min- 
imum, should be better compared with the EGRET 
subsample taken during the solar low level of activity. 
The neutron monitor rate [27] during the period of ob- 
servations by Fermi was, at least, 10% higher than the 
ground neutron rate values reported by EGRET Q at 
that time. We conclude therefore that there is a rea- 
sonable agreement with previous measurements if the 
solar activity is taken into account. 
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FIG. 4: Photon count map centered on Saturn position (marked with a circle) and coUected during the first 10 months 
starting from august 2008. The colored vertical scale is linear and smoothing has been applied to the image. The bin 
width used is 0.2°. The coordinates are celestial offset respect of Saturn position, the axes represent ecliptic coordinates. 



IV. OTHER SOLAR SYSTEM SOURCES 

A similar strategy has been used to search for the 
7-ray emission from the most interesting solar system 
objects, starting from the major planets. We have 
considered the photon count maps in the direction of 
major planets, including Jupiter and Saturn. For any 
selected object we have examined the count map to 
look for any evident excess of counts. No evidence 
of photon excess from the regions centered on these 
planets or within several degrees from these sources 
was found during the first 10 months of observations. 
Fig. m shows the count map of the events in celestial 
coordinates offsets relative to Saturn position. The 
search for the emission from populations of small solar 
system bodies. 15.] is currently in progress. 



lunar flux measured is higher than previous measure- 
ments but reasonably in agreement with the cosmic 
ray flux increase due the solar minimum activity in 
this solar cycle. We expect that with a larger sam- 
ple of events available after accumulating data for one 
year, Fermi should be able to explore the other fea- 
tures of the lunar spectrum, e.g. the tt*' peak and also 
to resolve the spatial structure of the emission from 
the lunar disk. We have started a search for the emis- 
sion from other bodies of the solar system objects too, 
in particular from Saturn and Jupiter regions. Dur- 
ing the first 10 months we don't have evidence of 7-ray 
emission in these regions. The search for the emission 
from populations of small solar system bodies is cur- 
rently in progress. 



V. CONCLUSIONS 

The gamma-ray emission from the Moon discov- 
ered by EGRET has been confirmed by Fermi and 
agrees in intensity for emission models that take into 
account the level of solar modulation. Our prelimi- 
nary flux estimation for the lunar 7-ray emission is 
F{E > lOOMey) = 1.1 ± 0.2 x 10"^ cm-^s-^ with 
a spectral index of —3.13 ± 0.03 obtained by fitting 
a simple power law between 100 MeV to 1 GeV. The 
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